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ABSTRACT 

We present an observational search for the possible correlation between cos- 
mic microwave background (CMB) polarization map and soft X-ray background 
(SXRB) based on the ROSAT All-sky Survey (RASS) archive and WMAP five- 
year observations. This is motivated by the fact that some of the CMB polar- 
ization may arise from the scattering of CMB photons due to the free electrons 
generated by X-ray heating sources in the epoch of reionization. Detection of 
such a correlation allows one to study the role of X-ray heating in the process 
of reionization. However, the cross angular power spectrum of the CMB polar- 
ization and SXRB maps constructed from ROSAT RASS and WMAP five-year 
maps is consistent with no correlation. We attribute this negative detection to 
both the extremely weak signals and large instrumental noises. While a future 
search is needed with high sensitivity instruments for both CMB polarization and 
soft X-ray photons, our current results can still be used as a useful constraint on 
the effect of X-ray heating in the epoch of reionization. 

Subject headings: cosmology: cosmic microwave background — large-scale struc- 
ture — X-rays: galaxies 



1. Introduction 



Many efforts have been made so far to search for the correlation between cosmic mi- 
crowave background (CMB) and X-ray sky, in an attempt to detect the signatures of 
both the integrated Sachs- Wolfe effect resulting from time-evolving potentials of large-scale 
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structures at late times, which can be traced by various X-ray sources, and the Sunyaev- 
Zel'dovich (SZ) effect due to the scattering of CMB photons by the ionized electrons in 
clusters of galaxies. Yet, no significant detections of the correlation signals h ave been re- 
ported using the ROSAT RASS data (IKneissl et al.l 1 199 71 ; iDiego et al.l 120031 ) . in contrast 
to the positive results reve aled by tracers of large-scale structures, including the hard (2 



10 keV) X-ray background (IBoughn fc Jahodalll993l ; iBennet et al.lll993l; iBandav et al. 
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source correlation generated at late times, CMB polarization map may also correlate with 
the unresolved soft X-ray background (SXRB) presumably originated from supernovae, X- 
ray binaries and quasars at high-redshift beyond z = 6. These sources are believed to 
play a signifi cant role in the heating of neutral gas through pho toionization at the epoch of 
reionization (IDijkstra et al.ll2004l ; IChen fc Mir alda-Escudel 120041 ). while the ionized electrons 
generate the CMB polarization signa ls, depending on b oth electron scattering optical depth 
and primordial quadrupole of CMB. iDore et al.l (120071) note tha t such polarization signals 
on small scales are on order of 0.01 fiK (see also lLiu et al.l 1200 ll . for theoretical prediction) 
. But as of today a quantitative estimate of the correlation signals is still unavailable in the 
literature, and any of such detections will challenge the existing instrumental limits. The 
goal of this study is to make the first attempt at searching for such a correlation, using the 
current upper limits of the resolved diffuse SXRB and the recent WMAP polarization map. 



A considerably large fraction (~ 94%, see iMoretti et al.l 120031 ) of soft and hard X-ray 
backgrounds have been resolved into discrete sources, and the maximum admitted ranges o f 
the unresolved flux at different energy bands have been summarized in IWu fc Xud (120011 ). 
If internal physical processes are properly treated, diffuse emission of hot gas in groups 
and clusters should also account for the residual SXRB after the removal off the discrete 



sources ( jXue fc Wul 120031 ). This indeed restricts the room for the presence of significant X 



ray e mission from high-redshift sources at the epoch of reionization (but see ISalvaterra et al. 



20071 ). So the fraction of soft X-ray which is emitted from high-redshift sources at the epoch 
of reionization could be not more than 5%. Furthermore, even if the correlation is detected, 
the signature may have been interlaced with those caused by the SZ effect, in the sense that 
the CMB polarization anisotropics may al so be generated by hot electrons in clusters and 
even groups of galaxies (jShimon et al.ll2006l ). unless foreground clusters and groups including 
distant ones can be masked in the search. In particular, the fluctuations of SXRB in the 
epoch of reionziation may appear at smaller scales comparable to or smaller than the sizes of 
clusters, because X-ray photons usually have relatively long free paths. This adds to further 
difficulties of our search for the signatures of X-ray reionization. Keeping these disadvantages 
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in mind, we perform an analysis of the correlation between the WMAP five-year polarization 
map and ROSAT RASS archive at energy band 0.2-2 keV, in the hope that current limits 
are helpful for placing robust constraints on models of X-ray reionization and for guiding 
future works on this topic. 



Data 



For WMAP five-year maps, we choose to work with the V and W bands in order to 
reduce foreground contamination. Furthermore, we restrict the analysis to the differently 
assembly polarization data denoted by VI, V2, Wl, W2, W3, and W4, which are called WP 
sky maps hereafter. Full five-year sky maps are obtained by performing a weighted, pixel- 
by-pixel, mean of the single-year maps. Polarization signals are characterized by the Stokes 
parameter Q, U, and V. Because Thomson scattering can not generate circular polarization, 
we have V = 0. WMAP observes the sky with two orthogonal l inear polarization m odes per 
feed, which are sensitive to the Stokes parameter I, Q, and U (IHinshaw et al.l 120071 ). These 
data were released in HEALPbfl format with different resolutions. The map used in this 
work is of Res 9 ( iV S id e = 512, which means the angular resolution is 6.87' ). 



The ROSAT RASS (jSnowden et al.l 1 199 71 ) covers the soft X-ray energy band ranging 
from 0.1 to 2 keV. Th e X-ray sky in lower energies is severely contaminated by local emis- 
sion (jSliwa et al.ll200ll ) while the X-ray sky in high energies shows a significant contribution 
from extragalactic AGNs. The RASS intensity (Jsxrb) and noise (a) maps are constructed 
by Jsxrb = (C — B)/t exp and a = a/C + B/t ey[p , respectively, where C is the photon counts, 
t exp is the exposure time and B denotes the contamination, which comprises sources, particle 
background, scattered solar X-rays, long-term enhancements, and so on. 

We use two types of SXR B maps constructed fr om different algorithms. The first one 
(hereafter RSI) is presented by lSnowden et al.l (119971 ). RSI covers approximately 98% of the 
sky in 1/4 keV, 3/4 keV, and 1.5 keV bands, each of which is further divided into two bands 
(R1+R2, R4+R5, R6+R7, respectively) with an angular resolution of about 12'. These maps 
were constructed using an equal-area polar projection and divided into six areas according to 
the six different projection pivots, in which the contamination sources, such as the particle 
background , scattered solar X-rays, and long-term enhancement, are excluded. 

To produce an SXRB map with the angular resolution compatible to that of WMAP 
data, we construct the SXRB map (hereafter RSII map) in HEALPix format with Res 9 
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(iV S ide = 512) from RASS archive with energy band ranging from 0.2 to 2 keV, in which each 
pixel covers an area of 6.87' x 6.87' in the sky. Except for the correction to exposure time, 
we make no attempt to remove the bright point and diffuse sources from the local universe 
in order to optimize the detection of cross-correlation with the WMAP polarization maps, 
even if it has a local origin. 



3. Statistical Method and Systematic Effects 



Following iKogut et al.l (120071 ). we expand the polarization fluctuations into generalized 
spherical harmonics 



pW(n) = Q(n)±iU(n) = ^a TA „ ^Y tm (n). 



The polarization fluctuations are conventionally decomposed into E- and B-modes. Intro- 
ducing a weighting function u> p (n), we can write the E-mode as 



i J dnw p (n)[Q(n)( 2 Y; m (n) +_ 2 Y/ m (n)) 

-ic/(n)( 2 y;jn) -_ 2 r; m (n))], 
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and its angular power spectrum is given by 
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The fluctuation of SXRB intensity AJ can be decomposed in spherical harmonics as 
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A/(n) = ^ a ^m^™( n ) 
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and the corresponding angular power spectrum is 
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We define the cross-correlation power spectrum between CMB E-mode map and SXRB as 



2£+l ' J m m ' 
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To construct the angular power spectrum from the observations, masks are often em- 
ployed to account for the incomplete or inhomogeneous sky coverages, or excisions of con- 
taminated regions, which are all incorporated in a window function W(n). For the WMAP 
polarization map, w e choose the window function W p (n) from the WMAP prescription 
([Bennett et al.l 120031 ). and we apply the natural weighting assumption to the window func- 
tion for the ROSAT map, namely Wi(n) = in the excised regions and Wi(n) = 1 elsewhere. 
The cross-power spectrum between SXRB and WMAP polarization maps Cf in the presence 
of incomplete sky coverage is modified by 
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Such a pseudo cross-power spectrum Cf can be related to the unbiased whol e sky power 
spectrum through their ensemble averages in terms of the algorithm MAS TER (IHivon et al. 
20021 ; lHansen fc G6rski|[2003l : Ichallinor fc ChorJl2~005l ; iBrown et al.ll2005h . 



where Mw describes the mode-mode coupling 
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and W is the cross-power spectrum of the window functions 



(9) 



(10) 



with w\ m and w[ m defined as the spherical harmonic coefficients. 

The angular power spectrum Cf should be further corrected for the smoothing effects 
of both instrumental beams {B), finite pixels (p) and noise (N) 



where B\, and Bf, are the Fourier transforms of beam patterns of ROSAT and WMAP, 
respectively. The pixel window function in Eq. (11) is defined by 
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and the m-averaged window function reads 

2£ + 



in which pe m is the spherical harmonic transform of pixel function P p i x (n) = l/fi P i x , and Q p i x 
is the solid angle subtended by one pixel. 

Finally, the angular power spectrum should be properly binned to reduce the system- 
atic effects described above, and the associated errors should also be estimated. With the 
algorithms MASTER and its extension to polarization, we can convert (Cf) to an 'unbiased' 
whole sky angular power spectrum C and evaluate the errors of C. In order to calculate the 
errors and covariance of the C, we produce a number of mock observational maps by Monte 
Carlo simulations. So the elements of the covariance matrixes are given by 

Cbb' = \ (p*> ~ (^b)Mc) \Cy — (Cb^McJy • (1^) 
The error bars on C\, are then given by the square root of the diagonal elements of C 

AC b = Cli 2 . (15) 



4. Results 

Before we proceed to the study of the cross-correlation between the CMB polarization 
map and SXRB, we present an extensive analysis of the auto-correlation angular power 
spectrum of the two backgrounds revealed by ROSAT and WMAP, respectively. This allows 
us to have a better understanding of the statistical properties of our data sets. 

We firstly work with the RSI maps in all the seven channels. Pixelizing the RSI sky under 
the angular resolution of 12' gives rise to a total pixels of iV S jde = 256 in the HEALpix format. 
In order to match the WMAP angular resolution, the data are somewhat oversampled so that 
a value of iV S ide = 512 is actually used, which has no influence on the construction of power 
spectra at large angular scales. We have tested a few of weighting functions to compensate 
for the incomplete sky coverage, and found that the results suffer from the contamination 
of foreground structures, such as the Galaxy and bright point sources in the soft X-ray 
band. We eventually use the R6 data within b > 40° and I G [70,250], the so-called MSR 
mask, which is suggested as a better repr esentative of the diffuse cosmological SXRB and 



less affected by the foreground structures (ISliwa et al.ll200ll ). The angular power spectrum 



is constructed using the HEALPix package and the result is shown in Figure 1(a). 
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Next, we apply the same technique to the RSII map with pixel resolution = 512. 
Note that this sample has a wider energy coverage ranging from 0.2 to 2 keV. Since we make 
no attempt to edit the raw data except correct for the exposure time, the RSII map looks 
very noisy and a large fraction of the sky is apparently dominated by the emission of the 
Galaxy. We demonstrate the resulting angular power spectrum with mask MSR in Figure 
1(b). 

Despite the fact that the two samples have different energy coverages , the overall angular 
power spectra revealed by RSI and RSII maps are similar in shape. Two samples differ by 
several orders of magnitude in amplitude because the RSI maps divide up the flux of SXRB 
among seven energy bands ,and only the power spectrum of R6 band map is plotted in the 
Figure 1(a). The slight excess in the RSII power spectrum at large scales below I ~ 100 can 
be attributed to the emission from the galactic plane of the Galaxy. Furthermore, the two 
angular power spectra both have a steeper slope than £{£ + l)Ce oc £ 0,8 , indicating that the 
SXRBs are actually dominated by local sources (i.e., the Galaxy) rather than X-ray galaxies 
distributed at cosmological distances. Recall that the clustering of galaxies contributes an 
angular power spectrum of roug hly £{l + l)C e oc £°- 8 , provid ed that the sources follow a 



power-law angular power spectrum with an index of ~ 0.8 (IScott fc White! Il999l ). Yet, 
such a local component makes no contribution to the cross-correlation of SXRB and CMB 
polarization to be explored below. 

Now we turn to the CMB polarization sky manifested by the Stokes Q and U parameters 
in the WP map. In a similar way to the analysis of the SXRB maps above, we apply two 
masks in the analysis of CMB polarization map; one is the full-sky WP data without an y 



masks, and another is the so-called P06 mask used by the WMAP team (IPage et al.l 120071 ). 
For the latter, 25.7% of the sky is masked near the Galactic plane and around other strong 
sources of contamination, eliminating the effect of synchrotron polarization by the Galaxy, 
dust, point sources, etc. Indeed, the angular power spectra of polarization map with and 
without the P06 mask look rather different from each other, as illustrated in Figure 2 using 
the CMB TE cross-power spectrum. Apparently the structured power spectrum of the full- 
sky map originates from the foreground, and the structures vanish immediately when the 
contaminated regions of the sky are masked with the P06 mask. Shown also in the Figure 
2 is the theoretically predicted CMB TE power spectrum calculated in terms of publicly 
available code CMBFASt|. It appears that the observationally measured and theoretically 
predicted power spectra agree each other nicely when the P06 mask is applied. 

Finally we carry out the study of cross-correlation between the CMB polarization sky 
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and SXRB using the ROSAT RSI and RSII samples and WP maps. Construction of cross- 
power spectra can be made straightforwardly, and the error bars of the cross-power spectra 
are estimated from 100 mock WMAP observational polarization maps and 100 mock SXRB 
maps using Eq. (|14p and Eq. (IT51) . The results of the cross-power spectra are shown in 
Figures 3 and 4 for RSI-WP and RSII-WP cross-correlations, respectively, for which we have 
considered the sky maps both with and without masks in order to demonstrate to what extent 
the angular power spectra are affected by the foreground sources. For the WP maps all the six 
polarization databases are employed and the results were averaged over the whole databases 
and the different observational band are illustrated, respectively. It turns out that, except for 
the RSII-WP cross-correlation without any masks, the rest are essentially consistent with the 
negative detection of cross-correlation between SXRB and CMB polarization maps. Recall 
that even in the full-sky coverage without the application of masks, such as MSR in the 
RSI (R6) sample, bright foreground X-ray sources have already been excised. Employment 
of the MSR allows one to further excise the Galactic disk from the samples. The absence 
of significant correlations between the RSI (with and without masks) or RSII (with mask 
MSR) samples and WP maps may therefore be of cosmological significance. We do detect 
a positive cross-correlation in the RSII-WP full-sky maps without any masks over angular 
scales of I ~ 200 — 400. Nevertheless, a careful examination of the RSII full-sky map shows 
that such a positive correlation is related to three strong X-ray emitting regions that are 
identified as Cyg Region, Gum Neb and Tau A in the WMAP point source catalog, which are 
strong X-ray source and supernova remnants. No correlation feature is detected once these 
three regions are excised or MSR and P06 masks are utilized (Figure 5). Actually, the result 
remains nearly unchanged if only one mask, either MSR or P06, is adopted to all the cases. 
The positive correlation associated with the three supernova remnants can be interpreted 
as the polarization of synchrotron in the regions or the polarization of CMB caused by the 
energetic electrons there, yielding the cross-correlation between X-ray emission and CMB 
polarization map. 



5. Discussion and conclusions 



The CMB polarization signals in large scale seen in the WMAP five-year data arise 
primarily from the scattering of free electrons in the epoch of reionization. At other scales, 
the primal quadrupole anisotropies are still dominant sources. But at small scales, the picture 
is more complicated. All types of the perturbations beyond linear order will contribute to 
the polarization signals. One possible source of polarization anisotropy at small s c ales is 
Thomson scattering by free electrons at the epoch of the reionization. iDore et al.l (12007) 
find such polarization signals may have rms of order ~ 0.01 fiK (see also iLiu et al.ll200ll ). 
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Other secondary contribution to the signals includes free electrons associated with large-scale 
structures through the warm-hot intergalactic medium (i.e. missing bary ons) in local universe 



(|Dave et al.ll200ll ). and very energetic electrons in clusters of galaxies (IBaumann &: Cooray 
20031 ) . These secondary effects (~ 10~ 4 — 10~ 5 fiK) are weaker than the polarization signals 
produced in the epoch of reionization and believed to be subdominant. Therefore, we focus 
on the possible reason for our negative detection of cross-correlation between ROSAT SXBR 
and WMAP CMB data. Indeed, there is good reason to believe that X-ray heating once 
playing a non-negligible contribution to the reionization of the universe. Primary sources of 
X-ray emission at high-redshift are high-mass X-ray binaries, supernovae and even QSOs. 
It is natural to expect a positive correlation between the CMB polarization and unresolved 
diffuse SXRB because both of them originate from the same epoch and follow the common 
underlying large-scale structures of the universe. 

The most likely explanation for the negative detection of the cross-correlation between 
CMB polarization map and SXRB obtained with the WMAP and ROSAT maps may be due 
to the extremely faint signals that are buried in the strong foreground and instrumentation 
noise. Indeed, both CMB polarization and X-ray emission from the epoch of reionization 
are rather weak, and the latter has not even been separated from current SXRB. But if 
we assume that the CMB pol arization signals due to the inh omogeneous process of the 



ionization of the universe (see Pore et al 



2007; Liu et al. 2001. for the estimation of such 



CMB polarization signals) were 100% correlated with the 5% of the soft X-ray emission 
which is emitted from high-redshift sources at the epoch of reionization, we can calculate 
this theoretical maximum correlations signals and take it as a rough estimate for the actual 
correlations signals. The results of the estimate, compared with the correlation results from 
the WMAP and ROSAT maps, are plotted as solid lines close to zero axes in Figure 3, 4, 
and 5. It is quit evident that even the theoretical maximum correlation signals are too weak 
to see in this measurement. Furthermore, if most of the CMB polarization is generated in 
the later time when the Lya heating, instead of the X-ray heating, dominates the process 
of reionization, the strength of the cross-correlation between CMB polarization and SXRB 
would be even weaker. Therefore, a quantitative estimate of the cross-correlation and its 
measurement errors with existing and future instruments will be desired to explain our 
negative detection ultimately, and to forecast further prospects. The prime difficulty in 
carrying out such a theoretical investigation is the uncertainty of the effect of X-ray heating 
sources and their cosmic evolution in the epoch of reionization. 

Obviously, such weak correlation signals will provide an ambitious target for upcom- 
ing experiments. Before a quantitative estimate is available, we can roughly evaluate the 
detectability of the correlation signals with the theoretical maximum signals estimate. Our 
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rough estimation of the noise cross-power spectra for Pl anck % (a = 2.3 fi K/5.0' fwhm, 0.85% 



of the sky) and ACT ( Atacama Cosmology Telescope; Kosowskyl bood ) or SPTPol (South 



Pole Telescope; iRuhl et all I2004T ) (a = 2.0 /j,K/I.T fwhm, 0.25% of the sky) polarization 
data correlating with ROSAT RASS data are illustrated in Figure 6, respectively. The levels 
of the noise cross-power for CMB experiments like ACT or SPTPol are above the theoretical 
maximum signals levels, while those for Planck are under the theoretical maximum signal 
levels. It appears that such theoretical maximum signals might be detected with Planck at 
£ < 900. 

In summary, if X-ray heating is one of the major sources to ionize the neutral hydrogen in 
the epoch of reionization, we may expect to find a positive correlation between the unresolved 
diffuse SXRB and CMB polarization maps. Intuitively, our negative detection based on the 
ROSAT and WMAP data indicates that such a correlation is rather weak, and experiments 
with higher sensitivities are required to make further explorations. Yet, our negative result 
can be used to constrain the theory of X-ray heating in the epoch of reionization. From our 
result, the theoretical prediction of the cross power spectra between R6 band SXRB and CMB 
polarization could not be out of range —0.8 < £(£ + l)/27r Ce < 0.7 (mK 10~ 6 cts/s/arcmin 2 ) 
at £ ~ 300, -1.5 < £{£ + l)/2ix C e < 1.2 at £ ~ 500, -3.3 < £{£ + l)/2vr C e < 3.2 at £ ~ 800 
(based on the result from RSI R6 band map and WP maps). And those between wider 
energy band (0.2 KeV ~ 2 KeV) SXRB and CMB polarization could not be out of range 
-3.3 < £(£+l)/2n C t < 3.2 (mK 10- 6 cts/s/arcmin 2 ) at £ ~ 300, -9.7 < £(£+l)/2n C e < 7.0 
at i ~ 500, -36.4 < £(£ + 1)/2tt C e < 30.2 at £ ~ 800 (based on the result from RSII map 
and WP maps). Meanwhile, a detailed study of the cross power spectrum would be is also 
worthy to guide future observations. 
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(a) (b) 




Fig. 1. — Angular power spectra revealed by the RSI and RSII samples. The la error bars 
are estimated using the WMAP simulations. Both angular power spectra are binned in terms 
of A£ = 30 
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Fig. 2. — TE angular power spectra constructed from WP maps(Vl, V2, Wl, W2, W3, W4 
data) with (b) and without (a) P06 mask, respectively. The solid lines represent the results 
obtained with the WP maps and the dashed lines are the theoretically predicted power 
spectrum generated by CMBFAST. The angular power spectrum constructed from the maps 
with masks are binned in terms of A£ = 50. 
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. Data: RSI R6 & WP maps 
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Fig. 3. — Cross-power spectra constructed from RSI R6 sample and WP six databases (dotted 
lines), and the open squares with error bars represent their averages over the whole database 
( top panel ), V ( middle panel ), and W ( bottom panel ) band. The solid line close to the 
zero axis represents the estimate of maximum correlation signals (see Section [5] for details). 
The small graph sketches the estimate correlations signals with magnifying Y-range. 
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Fig. 4. — Same as Figure 3, but for RSII sample. 
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Data: RSII & WP maps 
- Sky-coverage: full-sky & P06 mask '. 
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Data: RSII & WP V band maps , 
- Sky-coverage: full-sky & P06 mask 
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Fig. 5. — Same as Figure 3, but for the full-sky RSII sample with (a) the P06 mask on WP 
maps and (b) removal of three bright X-ray sources in Galactic plane (Cyg Region, Gum 
Neb and Tau A). 
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Fig. 6. — Theoretical maximum correlation signals one could expect compares with the es- 
timated noise of the cross-power spectrum for the upcoming experiments. The heavy solid 
lines represent the estimated maximum signals from CMB polarization map correlating with 
RSI R6 band map (a) and RSII map (b), respectively. The estimation of the noise cross- 
power spectrum for Planck polarization and SXRB maps is shown as dotted lines, while the 
estimation for the polarization observation like ACT or SPTPol and SXRB maps is shown 
as dash lines. 



